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Abstract 

Theoretical investigations on electrorheological (ER) fluids usually rely on 
computer simulations. An initial approach for these studies would be the 
point-dipole (PD) approximation, which is known to err considerably when 
the particles approach and finally touch due to many-body and multipolar 
interactions. Thus various work attempted to go beyond the PD model. Being 
beyond the PD model, previous attempts have been restricted to either local- 
field effects only or multipolar effects only, but not both. For instance, we 
recently proposed a dipole-induced-dipole (DID) model which is shown to 
be both more accurate than the PD model and easy to use. This work is 
necessary because the many-body (local-field) effect is included to put forth 
the many-body DID model. The results show that the multipolar interactions 
can indeed be dominant over the dipole interaction, while the local-field effect 
may yield an important correction. 
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I. INTRODUCTION 



An electrorheological (ER) fluid consists of polarizable particles in a nonconducting host 
fluid. If an external electric field is applied to an ER fluid, the particles aggregate and 
form chains parallel to the applied field. These chains may further aggregate to columnlike 
structures [|T|. To discuss ER effect, early theory is often based on the point-dipole (PD) 
approximation, e.g. H. Also, the PD approximation is often adopted in computer simulation 
H because it is simple and easy to use. Since many-body and multipolar interactions 
between particles have been neglected, the predicted strength of ER effect is of an order lower 
than the experimental results. Hence, much work has been done to sort out more accurate 
models @||H,[1. Although those methods are accurate, they are relatively complicated to 
use in dynamic simulation of ER fluids || . 

In a recent paper, we put forth the dip ole- induced dipole (DID) model to improve the 
point-dipole (PD) model |J. In that work, we considered pair interaction between polar- 
ized dielectric particles through the multiple image formula, but neglected the many-body 
interactions because it was believed that the multipolar interactions can be dominated over 
the many-body (local-field) effects. The claim was based on the calculation for a pair of 
particles. We should really question the quantitative accuracy of this claim, since in ER 
fluids chain and sheet structures are known to occur. In a long particle chain the local field 
at a dipole site is much larger than for a pair of dipoles, whereas for multipolar interactions 
this is probably not the case. In fact, this issue has already been quantified, for particles 
with large dielectric contrast with the liquid phase, by Martin and Anderson ||10|| . Thus we 
have to do some additional analysis to validate the claim. 

As an initial model, we may adopt an effective dipole factor in the multiple image formula 
||. More precisely, we regard that each polarized particle in an ER suspension is embedded 
in an effective medium with an effective dielectric constant e e . Thus, the usual dipole 
factor b — (ex — 62) /(e\ + 262) in the multiple image theory should be replaced by b' = 
(ei— e e )/(e!+2e e ) [[□]], where e e is the effective dielectric constant of the surrounding medium, 
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which is conveniently calculated in the Maxwell-Garnett approximation (MGA). Here e\ 
and e 2 are the dielectric constants of the suspended particles and host fluid, respectively. 
However, the situation is further complicated by the fact that the suspended particles in 
an ER fluid usually form anisotropic structures, e.g. chain and sheet structures, along 
the applied field and that e e can be anisotropic. Thus, e e should be calculated within the 



anisotropic MGA |12| , |13| . A preliminary calculation shows that e e should be increased and 
b' should be decreased as compared to b. Thus the DID interaction would be reduced. But 
of course the actual reduction should depend on the local-field factor of the anisotropic 
structures. Thus the many-body effects have been included through the volume-fraction 
dependent e e . 

So, in the presence of the surrounding particles, the interaction between two particles will 
be modified. By means of an effective dielectric constant, we can derive an analytic formula 
for the many-particle DID model. And more importantly, we can access the importance of 
the many-body (local-field) effects against the multipolar interactions. We believe this work 
is necessary because being beyond the fixed dipole model, the previous work dealed with 
either local-field effects only [|TJj] or multipolar effects only, but not both. 

II. FORMALISM 

We concentrate on the case where highly polarized dielectric particles of diameter d, 
dielectric constant e\ are embedded in a host fluid of e 2 . 

In the dilute limit, the dipole factor for an isolated particle may be given by 

6 = ii^. (i) 

On the other hand, based on the multiple image method || the dipole factor for a pair of 
particles with separation s has the following form: 

sinh a 3 
n L sinh(n + l)a 



&£. = 6 £ (-&)"[- 
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bi-bj:(2br[ -^— r, ( 2 ) 



n=0 



sinh(ra + l)a 



for transverse field (T) and longitudinal field (L) case, respectively, where a satisfies the 
relation cosh a = s/d. 

When a suspension containing many dielectric particles is subjected to an intense electric 
field, the induced dipole moments may cause the particles to form chains along the applied 
field, resulting in complex anisotropic structures. In this case, we may invoke the Maxwell- 



Garnett approximation (MGA) for anisotropic composites |12l , |T3| to obtain the effective 
dielectric function e e of the system. For the transverse field case when the electric field is 
applied perpendicular to the uniaxial anisotropic axis, the MGA has the form 

f^FP-, (3) 



PrteT + (3 - t3 T )e 2 ei + 2e 2 
whereas for a longitudinal field case when the field is applied along the uniaxial anisotropic 
axis, the MGA reads 

£eL ~ 62 = f 61 ~ 62 U) 

/W + (3 - f3 L )e 2 e 1 + 2e 2 ' y 1 

where (3t and (3l denote the local field factors perpendicular and parallel to the uniaxial 
anisotropic axis, and / the volume fraction of particles. These local field factors are defined 
as the ratio of the local field in the particles to the Lorentz cavity field |Tj|. For isotropic 



composites, /3l = /3r = 1, while both (3l and (3t will deviate from unity for an anisotropic 
distribution of particles in composites. These (3 factors have been evaluated in a tetragonal 



lattice of dipole moments [I2|] and in various field-structured composites [nj and they satisfy 
the sum rule f3i + 2/5^ = 3. In this case, the dipole factor for an individual particle may 
have this form [|ll]] 

61 " 6e (5) 



ei + 2e e ' 

Note, to obtain this equation, we have replaced e 2 in Eq.(|l|) with e e . That is, we assume 
the particle to be embedded in an effective medium. Thus the many-body (local-field) effect 
is included already in Eq.(|J). So far, to put forth the many-body DID model, we may 
consider a pair of dielectric spheres as well. The spheres are placed in an effective medium. 



A constant electric field Eo = Eqz is applied to the spheres, which contribute to each sphere 

a dipole moment given by pi and P2o(= Pio = £ e Eod 3 b' /8). Thus the multiple image formula 

can be developed in a similar way 0. The dipole moment pio induces an image dipole pn 

in sphere 2, while pn induces yet another image dipole in sphere 1. As a result, multiple 

images are formed. Similarly, p 2 o induces an image P21 inside sphere 1, and hence another 

infinite series of image dipoles are formed. Then we may obtain the sum of dipole moment 

inside each sphere, and hence the desired expressions for dipole factors: 

sinh a , 3 
/r () L sinh(n + l)a 



b%=b'j2(-b'n- 



6 " = %^n^^ 3 («) 

for transverse field and longitudinal field case, respectively. The two equations are nontrivial 
results indeed, which include the multiple image as well as many-body (local-field) effect. 

To discuss the interparticle force, we may take one step forward to calculate the force 
between two dielectric spheres. We have already obtained the dipole factors b'* and b* and 
hence the dipole moments, and the force can be calculated by an energy approach. In this 
case, the dipole energy (En) of a pair of particles may be determined by the dot product of 
dipole moment and the electric field, and hence the force between a pair of particles is given 
by the derivative of dipole energy with respect to separation, namely — dEn/ds. So, based 
on this relation, we may obtain the expressions for interparticle forces, respectively: 



with 



F* = -(l/8)e 2 E 2 d 3 bJ2(-b) n <S>, 

n=0 

oo 

F; = -(l/8)e#&E(2r$, 

n=0 

oo 

F't = -a/8)e eT E 2 d% £(-^) n $ , 

n=0 

oo 

F'l = -(l/8)e ei £ W L J2(2b' L r$, (7) 

n=0 



sinh a cosh a sinh[(ra + l)a] — (n + 1) sinh 2 a cosh[(n + l)a] , 
= (rf/3)sinh 4 [(n + l)«] ' ( } 
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where F£ (F£) indicates the interparticle force between a pair of particles in fluid host for 
transverse (longitudinal) field case, and F% (F'* L ) in an effective medium for transverse 
(longitudinal) field case. It is shown that the multipolar effect as well as many-body effects 
have been taken into account. Setting n up to 1 will yield directly the point-dipole forces, 
namely F T and F L . 

So far, we have put forth many-body DID model. We are now in a position to do some 
numerical calculations to discuss the (3 effect on reduction factor R = b'*/b*, which indicates 
the correction of the many-body (local-field) effect on the multiple image effect, as well as 
on the ratio of b* or b'* to b. Then, to investigate such effect on interparticle force, we also 
investigate the interparticle force normalized by point-dipole force. 

III. NUMERICAL RESULTS 

For numerical calculations, we choose e\ = 10e , e 2 = 2.5e , d = 2.0 xlO -7 and s/d = 1.01, 
where e is the dielectric constant of free space. From the numerical calculations, we find 
that a domination of the multipolar interactions is possibly over the many-body (local-field) 
effects. More precisely, in Fig.l, for the transverse field case, the correction due to the local- 
field effect is always small, because the reduction factor is close to unity. However, for the 
longitudinal field case, it is evident that the local-field effect can be large, especially at high 
volume fractions and/or large (3l- In which case, the correction due to local field cannot be 
neglected. In Fig. 2, we compute the ratio of b* or b'* to b versus the local field factor (3 L for 
different volume fractions. In Fig. 3, we discuss the dependence of (3l on the interparticle 
force normalized by point-dipole force for different volume fractions. From Figs. 2 and 3, it is 
obvious that the multipolar interactions can indeed be dominant over the dipole interaction, 
while the local-field effect may yield an important correction. 

The force expression is reasonable, as the n = term of Eq.(8) vanishes while the 
n = 1 term gives the correct dipole interaction force r~ 4 dependence. The divergence of the 
summation may be found for the longitudinal field case, due to the small distance s/d = 1.01 
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being used. 

From the results, it is clear that the correction due to many-body (local-field) effects can 
be very large, especially at large volume fractions. Nevertheless, for small volume fractions, 
the correction due to the local field can be neglected. The reduction of the magnitude of the 
effective interaction between two particles can be understood in a simple geometry in which 
the two particles interact in the presence of a third particle. An average over all possible 
positions of the third particle gives the desired reduction which is valid both for longitudinal 
and transverse fields. So, we believe our results are correct. 



IV. DISCUSSION AND CONCLUSION 

Here a few comments are in order. The ER effect arises from the polarization forces 
produced by a mismatch in the dielectric constants of the disperse and continuous phases. 
Typical ER fluids contain a disperse phase with volume fractions in the range of 0.05 ~ 0.50. 
In this sense, it is reasonable to set / = 0.1, 0.2, 0.3 in our calculations. On the other hand, 
we set €2 = 2.5eo, which is actual value for silicon oil and commonly employed as the 
host fluid. Also, the disperse phases are generally composed of solid, nonconducting, or 
semiconducting materials with dielectric constants in the range 2 ~ 40. Thus, the choice 
ei = 10eo is reasonable in our numerical calculations. In addition, the longitudinal local 
field factor (3l should satisfy 1 < (3l < 3, in order to ensure that ftr > 0. Regarding the 
separation between the pair of particles, various different values have been used as well, e.g. 
s/d = 1.03 etc., and similar results are obtained (not shown here). However, for the case 
where the separation is too large, e.g. s/d > 2, the multiple image effect will become small 



enough to be neglected [15 



Our many-body DID model may be used in various studies of the behavior of ER fluids, 
e.g., in a theory describing the ground state of ER solids or in a computer simulation of 
ER fluids. By including the many-body effects, the DID model can be used with a higher 
accuracy. Lastly, it is straightforward to extend our model to deal with polydisperse ER 
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fluids, and to discuss the magnetorheological suspensions. 

In summary, we have developed a many-body DID model by including the many-body 
(local-field) effect, to compute the interparticle force for an ER fluid. We applied the formal- 
ism to a pair of touching spherical particles embedded in an effective medium, and calculated 
the force as a function of the longitudinal local-field factor. The result show that the multi- 
polar interactions can indeed be dominant over the dipole interaction, while the local-field 
effect may yield an important correction, especially at larger volume fractions. 
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FIGURES 

FIG. 1. The reduction factor versus local field factor (3l for different volume fractions, / = 0.1, 
0.2, and 0.3. 

FIG. 2. The ratio of b* or b * to b versus local field factor 0l for different volume fractions, 
/ = 0.1, 0.2, and 0.3. 

FIG. 3. The ratio of interparticle force to point-dipole force for different volume fractions, 
/ = 0.1, 0.2, and 0.3. 
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